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Abstract 

It is of increasing contemporary interest to preserve ambient isotopy during 
geometric modeling. Bezier curves are pervasive in computer aided geometric 
design, as one of the fundamental computational representations for geomet- 
ric modeling. For Bezier curves, subdivision algorithms create control poly- 
gons as piecewise linear (PL) approximations that converge under Hausdorff 
distance. A natural question is whether subdivision produces topologically 
reliable PL approximations. Here we focus upon ambient isotopy and prove 
that sufficiently many subdivisions produce a control polygon ambient iso- 
topic to the Bezier curve. We also derive closed-form formulas to compute 
the number of subdivision iterations to ensure ambient isotopic equivalence 
in the resulting approximation. This work relies upon explicitly construct- 
ing homeomorphism and ambient isotopy, which provides more algorithmic 
efficiency than only showing the existence of these equivalence relations. 
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1. Introduction 



Ambient isotopy is a stronger notion of equivalence than homeomorphism. 
Precisely, isotopy is a homotopy with the requirement of a homeomorphism 
at each value of the time parameter, which is particularly applicable for 
time varying models. Based on the concurrent result of angular convergence 
and homeomorphism [15] , here we establish the stronger equivalence relation 
between Bezier curves and the control polygons. 



A Bezier curve (Definition 3.1) is characterized by an indexed set of 



points, which form a PL approximation of the curve, called a control polygon 



(Definition 3.1 ). The de Casteljau algorithm [HldS] is a subdivision algorithm 
associated to Bezier curves which recursively generate control polygons more 
closely approximating the curve under Hausdorff distance [2"3j . 

There is contemporary interest [TJ [2j [18j [22] to preserve topological char- 
acteristics such as homeomorphism and ambient isotopy between an initial 
geometric model and its approximation. In particular, a control polygon am- 
bient isotopic to the Bezier curve is considered having the same embedding in 
Knot theory. These are fundamental for applications in computer graphics, 
computer animation and scientific visualization. 

However, there may be substantial topological differences between Bezier 
curves and their control polygons. There is an example showing an unknotted 
Bezier curve with a knotted control polygon [3j ETJ . An example of a knotted 
Bezier curve with an unknotted control polygon was constructed recently [H] . 

Although there may be substantial topological differences between a Bezier 
curve and its initial control polygon, we prove that sufficiently many subdivi- 
sions will generate a control polygon that is ambient isotopic (and therefore, 
necessarily homeomorphic) to the curve. We provide closed-form formulas to 
compute a sufficient number of subdivisions to achieve this ambient isotopy. 

2. Related Work 

The topological characteristics for a Bezier curve and the control poly- 
gon have been studied in the literature. There is an existing proof [2S] of 
homeomorphism under subdivision relied upon the hodographP] and did not 
provide the number of subdivision iterations. A geometrically constructive 



2 The derivative of a Bezier curve is also expressed as a Bezier curve, known as the 
hodograph [§]. 
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proof is established simultaneously, which provides the number of subdivision 
iterations |15j . Topologically reliable approximation in terms of homeomor- 
phism of composite Bezier curves was also studied by algorithmic techniques 
that do not completely rely upon de Casteljau algorithm. Instead, that work 
starts with choosing so called "significant points" [BJ- 

Ambient isotopy for geometric models was obtained by a certain algorithm 
distinct from subdivision [TT]. Recent papers pEJ [16] present algorithms to 
compute the isotopic polygonal approximation for 2D algebraic curves. Com- 
putational techniques for establishing isotopy and homotopy have been estab- 
lished regarding algorithms for point-cloud data by "distance-like functions" 
[5]. The input data here is different of a Bezier curve, with the important 
additional result of the number of subdivision iterations needed to attain an 
ambient isotopic approximation of the original Bezier curve. 

Ambient isotopy under subdivision was previously established (22] for 
3D Bezier curves of low degree (less than 4), where a crucial unknotting 
condition was trivial. The results presented here extend the earlier theorem 
to Bezier curves of arbitrary degrees. The focus on higher degree versions 
was motivated by applications in molecular simulation where Bezier curve 
models are created on input of hundreds of thousands of points, with interest 
in having curves that are at least C l . Preserving that continuity over low 
degree models on this magnitude of points would be extremely tedious. 

Recent progress regarding isotopy under certain convergence criteria has 
been made pD, [13]. These are nice theorems applicable to Bezier curves. 
However, they do not necessarily provide formulas to compute the number 
of subdivision iterations, which is computationally crucial, as Bezier curves 
are used in many practical areas. 

Convergence of the control polygon to a Bezier curve under subdivision 
in terms of Hausdorff distance is well known (2U [25]. For a space curve, it is 
known how to construct a "nice" tubular neighborhood, that is, a neighbor- 
hood bounded by a nonsingular pipe surface [IT]. Based on the convergence 
in distance, one may obtain a control polygon that fits inside the tubular 
neighborhood. Pipe surfaces have been studied since the 19th century |20j, 
but the presentation here follows a contemporary source p2]. These con- 
temporary authors [17J provide sufficient conditions for a pipe surface to be 
nonsingular when its spine curve is regular and C l , as is true for B. These 
authors perform a thorough analysis and description of the end conditions 
of open curves. The junction points of B are merely a special case of that 
analysis. 
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3. Definitions and Notation 

Mathematical definitions, notation and a fundamental supportive theo- 
rem are presented in this section. More specialized definitions will follow in 
appropriate sections. The standard Euclidean norm will be denoted by || ||. 



Definition 3.1. The parameterized Bezier curve, denoted as B(t), of de- 
gree n with control points P m G M 3 is defined by 



Bit) = J2 B m,n(t)Pm,t G [0,1], 



m=0 



n 



where B m ^ n {t) = [J t m {l — t) n 111 and the PL curve given by the points 

{Pq, P\, . . . , P n } is called its control polygon. When Pq = P n , the control 
polygon is closed. Otherwise when P ^ P n , it is open. 

In order to avoid technical considerations and to simplify the exposition, 
the class of Bezier curves considered will be restricted to those where the 
first derivative never vanishes. 

Definition 3.2. A differentiate curve is said to be regular if its first deriva- 
tive never vanishes. 

Definition 3.3. A curve is said to be simple if it is non- self -intersecting. 



The Bezier curve of Definition 3.1 is typically called a single segment 



Bezier curve, while a composite Bezier curve is created by joining two or 
more single segment Bezier curves at their common end points. 

We use B to denote a simple, regular, C 1 , composite Bezier curve 
in IR 3 , throughout the paper. 



Definition 3.4. [21] Let X and Y be two non-empty subsets of a metric 
space (M,d). We define their Hausdorff distance fi(X,Y) by 

u(X,Y) := maxjsup inf d(x, y),sup inf d(x,y)}. 
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Figure 1: An exterior angle 



Exterior angles were denned [19] in the context of closed PL curves, but 
are adapted here for both closed and open PL curves. Exterior angles unify 
the concept of total curvature for curves that are PL or differentiable. 

Definition 3.5. The exterior angle between two oriented line segments, 
denoted as P m _\P m and P m P m+ i, is the angle between the extension oj ' P m _\P m 



and P m P m+ i, as shown in Figure 1(a), Let the measure of the exterior angle 
to be a m satisfying: 

< a m < 7r. 

This definition naturally generalizes to any two vectors, V\ and Vi, by joining 
these vectors at their initial points, while denoting the measure between them 
as T](vi,V2), as indicated in Figure 2(b). 

Definition 3.6. The curvature of a C 2 curve C(t) parametrized on [a,b] 
is given by 

\\C'(t) x C"(t)\\ 

«(*) = " \c'(tw tG M- (i) 

Its total curvature [£j is the integral: dt. 

Total curvature can be defined for both C 2 and PL curves. In both cases, 
the total curvature is denoted by T K (-). The unified terminology is invoked in 
Fenchel's theorem, which follows and is fundamental to the work presented 
here. 

Definition 3.7. [TS] The total curvature of a PL curve in IR 3 is the sum 

of the measures of the exterior angles. 
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Fenchel's Theorem [5] presented below is applicable both to PL curves 
and to differentiable curves. 



Theorem 3.1. [U Fenchel's Theorem] The total curvature of any closed 
curve is at least 2k, with equality holding if and only if the curve is con- 
vex. 

Denote a PL curve with vertices {Pq, Pi, ■ ■ ■ , P n } by P, and the uniform 
parametrization [23] of P over [0, 1] by Z(P)[o,i]. That is: 

l(P) m ( 3 -) = Pj for j = 0,1,- ■■ ,n 
and Z(P)[o,i] interpolates linearly between vertices. 

Definition 3.8. Discrete derivatives [23] are first defined at the param- 
eters tj = -, where 

l(P) [0A] (t J )=P 3 

for j — 0, 1, • • • , n — 1. Let 

p; = np) m ^) = ^^. 

Denote P' = (Pg, P{, • • • , P' n _i)- Then define the discrete derivative for '/(P)[o,ii 
as: 

l'(P) m =l(P') m . 

Intuitively, the first discrete derivatives are similar to the slopes defined 
for univariate real- valued functions within an introductory calculus course. 

The definition given here is specialized to M. n as sufficient for the appli- 
cations intended, but further generalizations are available [TO] : 



Definition 3.9. Two subspaces ofW n , denoted by X and Y, are said to be 
ambient isotopic if there exists a continuous map H : W 1 x [0, 1] — > R n 

satisfying: 

(1) H(-,0) is the identity, 

(2) H(X, 1) = Y, and 

(3) Vt G [0, 1], H(-,t) : lR n lR n is a homeomorphism. 
The map H is an ambient isotopy between X and Y . 



6 



Definition 3.10. The pipe surface of radius r of a parameterized curve 
c(t), where t G [0, 1] is given by 

pit, 6) = c(t) + r[cos(9) n{t) + sin{6) b(t)], 

where 9 e [0, 2tt] and n(t) and b(t) are, respectively, the normal and bi- 
normal vectors at the point c(t), as given by the Frenet-Serret trihedron. 
The curve c is called a spine curve. 

Remark 3.1. The paper J77| / provides the computation of the radius r only 
for rational spline curves. However, the method of computing r is similar for 
other compact, regular, C 2 , and simple curves, that is, taking the minimum 
ofl/n max , d m i n , andr en d, where K max is the maximum of the curvatures, d m i n 
is the minimum separation distance, and r en d is the maximal radius around 
the end points that does not yield self-intersections. 

4. Construction of Homeomorphisms 

Homeomorphism is a necessary condition for ambient isotopy. The exis- 
tence of a homeomorphism between E and the control polygon is shown in 
the companion paper [15] . However, constructing the ambient isotopy here 
relies upon explicitly constructing a homeomorphism. The explicit construc- 
tion provides more algorithmic efficiency than only showing the existence of 
these equivalence relation^} 

We shall find two criteria {Hypotheses 1, 2) to ensure the ambient isotopy 
between a PL curve and a C 2 curve. We then show that these criteria 
can be achieved for a Bezier curve by subdivision. Based on this, we shall 
construct an ambient isotopy between a Bezier curve and the control polygon 
in Section [U 

Let L be an open, compact, and simple PL parametric curve and C to 
be an open, compact, simple, regular, and C 2 polynomial curve of degree n 
in IR 3 . Let L(t) : [0, 1] -»■ 1R 3 and C(t) : [0, 1] -»■ 1R 3 be the corresponding 



3 Thcre is an alternative proof of the ambient isotopy |12j . which requires the convex 
hulls of sub-control polygons to be contained in the tubular neighborhood determined 
by a pipe surface. This containment makes the proof easier, but likely needs too many 
subdivision iterations and produces too many PL segments. However it is valuable for 
theoretical insight. 
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parameterizations of L and C. We consider the homeomorphism between L 
and C, when L(0) = C(0), L(l) = C(l), that is, LUC forms a closed curve^J 
Since C is a compact, simple, regular, and C 2 curve, there exists a non- 



singular pipe surface (Remark 3.1) of radius r. Denote the disc of radius 
r centered at C(t) and normal to C as D r (t). Let a pipe section to be 
T = U t6 r 1 j D r (t). Denote the interior as int(T), and the boundary as dT. 
Note that the boundary dT consists of the nonsingular pipe surface (Defini- 
tion 3.10) and the end discs D T (0) and D r (l). 

Definition 4.1. Define 6{t) : [0, 1] -> [0, tt] 6?/ 

e(t) = r,{C(t),L'(t)), 



where the function T](-, •) was previously given in Definition 3.5 



The following Hypotheses 1 and 2 will play a central role to 
the following lemmas, and in establishing the homeomorphism. 

Hypothesis 1: L \ {L(0),L(1)} C int(r); and 
Hypothesis 2: T K (L) + max^^i] 6{t) < |, 

where V is the pipe section of C and T K (L) denotes the total curvature of L. 

Lemma 4.1. Suppose L is a sub-control polygon and C is the corresponding 
Bezier sub-curve. Then Hypotheses 1 and 2 can be achieved by subdivision. 

Proof: By the convergence in Hausdorff distance under subdivision, suf- 
ficiently many subdivision iterations will produce a control polygon that fits 
inside a nonsingular pipe surface. Furthermore, by the Angular Convergence 
[T5l Theorem 4.1] and the lemma [151 Lemma 5.2] in the companion paper, 
possibly more subdivisions will ensure that each sub-control polygon lies in 
the corresponding nonsingular pipe section, which is the Hypothesis 1. De- 
note the number of subdivision iterations to achieve this by i\. 

By the Angular Convergence [15, Theorem 4.1], T R (V) converges to 
under subdivision. Because the discrete derivative of the control polygon 
converges to the derivative of the Bezier curve [23] under subdivision, 8(t) 



'Any pairwise disjoint simple open curves are trivially ambient isotopic. 
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converges to for each t G [0,1]. So Hypothesis 2 will be achieved by 
sufficiently many subdivision iterations, say i<i- (The Details to find i\ and 
6 2 are in Section [6|) □ 



Remark 4.1. When we consider the number of subdivision iterations later, 
we shall show that the sufficient number of subdivisions for Hypothesis 2 is 



at most one more than that for a weaker condition T K (L) < |. (Remark 6.1 
in Section^ 

Remark 4.2. To obtain some intuition for these hypotheses, restrict our at- 
tention to a Bezier curve. Consider L to be a sub-control polygon and C to 
be the corresponding sub-curve. Hypothesis 1 will ensure that L lies inside a 
nonsingular pipe section, while Hypothesis 2 will ensure a local homeomor- 
phism between L and C . In particular, Hypotheses 1 and 2 will be sufficient 
for us to establish the one-to-one correspondence using normal discs of C . 

Hypotheses 1 and 2 are assumed in the rest of the section. 

Define a function L(t) : [0, 1] — > L by letting 

L(t) = D r (t) n L, (2) 

where D r (t) is the normal disc of C at t. 

Define a map h : C — >■ L for each p e C by setting 

h(p) = L{C-\p)). (3) 

We shall show that h is a homeomorphism. The subtlety here is to demon- 
strate the one-to-one correspondence by showing each normal disc of C in- 
tersects L at a single point (which will be the main goal of the following), 
and intersects C at a single point (which will be easy), under the assumption 
of Hypotheses 1 and 2. 

4-1. Outline of the Proof 

For an arbitrary to £ [0, 1], the associated normal disc is denoted as 
D r (to). Following is the sketch of proving that D r (to) intersects L at a single 
point. (See Figure 12} ) 



1. The essential initial steps are to select a non-vertex point of L, denoted 
as w, a plane, denoted as Q, and an angle, denoted as 9(to): 
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Slope of h = L'(to) 




QHD»(to) n fi HC(to) 



Figure 2: Each normal disc intersects L at a single point 



(a) Define w and Q: Pick a line segment of L whose slope is equal 
to L'(t ), denoted as h. Choose an interior point of h, denoted as 
w. Let Q be the plane that contains w and is parallel to D r (t ). 
(We use w to define two sub-curves of L, a 'left' sub-curve which 
terminates at w, denoted as Li, and a 'right' sub-curve which 
begins at w, denoted as L r .) 

(b) Consider rj(C (t ) , L' (t )) = 6(t ) (Definition gj). Since Q is par- 
allel to D r (t ), a normal vector of Q, denoted by % has the same 
direction as C'(t ) and r)(nn,fo) = r)(C'(to), L'(t )) = 0(t ). 




Figure 3: Similar angles 9(t 



Remark: Since rj(nn, h) = 9(to), Hypothesis 2 implies that T K (L)+T](nn, h) = 
T K (L)+Q(to) < |. Since w is an interior point of h, the angle determined by Q 
and L/, and the angle determined by and L r , have the same measure 0(to), 
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as shown in Figure[3| So we obtain the similar inequalities T K (Li)+0(to) < 



and T K (L r ) + 6 (to) < | , which will be crucial. 

Prove, by Hypothesis 2, that QC\L r = w . Similarly, show that flHLi 
w. So Q fl L = w. (Lemma 4.3) 



3. Prove that any plane parallel to f2 intersects L at no more than a single 



point. (Lemma 4.4) 



Since D r (to) || fi, it will follow that D r (to) intersects L no more than a 
single point. Show, using Hypothesis 1, that D r (to) must intersect L, 



and hence D r (to) fl L is a single point. (Lemma 4.5) 



Preliminary Lemmas for Homeomorphisms 
In order to work with total curvatures of PL curves, an extension of the 



spherical triangle inequality [27], given in Lemma 4.2, will be useful, similar 
to previous usage by Milnor |19j . 




B 



Figure 4: Spherical triangle AABC 

Spherical triangle inequalities: Consider Figure [4| and the three 
angles ZAOB, ZBOC, and /LAOC ', formed by three unit vectors OA, OB, 
and 0(5. (Note the common end point O. When we consider angles between 
vectors that do not share such a common end point, we move the vectors to 
form a common end point.) Denote the arc length of the curve from A to B 
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as £(AB), and similarly for that from B to C as £(BC) and that from A to C 
as £(AC). The triangle inequality, £(AB) < £(BC)+£(AC), of the spherical 
triangle AABC provides that 

ZAOB < ZBOC + ZAOC. (4) 

Lemma 4.2. Suppose that v\, v 2 , ■ ■ ■ > ^m; where m G {3, 4, . . are nonzero 
vectors, then 

r)(vi,v m ) < T](vi, v 2 ) + r)(v 2 , v 3 ), + +r)(v m -i, v m ). (5) 

Proof: The proof follows easily from Inequality |4| □ 
Now, we adopt the notation shown in Figure [2] and formalize the proof 
outlined in Section [47L We assume that the sub-curve on the right hand side 
of Q in Figure [3] is L r , and the other one is Li, where we denote the set of 
ordered vertices of L r as 

{v ,v u . . .,«„}, 

with Vq — w. 

By Definition 4.1, we have 6(t ) < max te [ ,i] 0(t). It is trivially true that 
T K (L r ) < T K (L), so that with Hypothesis 2: T K (L) + max te [o,i] 9(t) < |, we 
have 

T K (L r ) + 9(t Q ) < T K {L) + max 9(t) < |. (6) 



The statement and proof of Lemma 4.3 depend upon the point w chosen 
in Step 1 of the Outline presented in Section 4.1 There, the point w was 
defined as an interior point of a line segment h~ of L, so that w is precluded 
from being a vertex of the original PL curve L. 

Lemma 4.3. The plane Q intersects L only at the single point w. 

Proof: Here we prove Qr\L r = w. A similar argument will show QdLi = w. 
The oriented initial line segment of L r is wv\ which lies on h. So 

r}(n Q , wvi) = r](n n , h) = 9(t ) < -. 

For a proof by contradiction, assume that Q intersects L r at some point 
u other than w. The possibility that wv\ C Q is precluded by 9 (to) < 7r/2, 
so the plane Q intersects wv\ only at w. So u £ wv\. 
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Figure 5: The intersection u generates a closed PL curve 

Denote the sub-curve of L r from w to u as L(wu). Then, since u ^ wv\, 
the union, L(wu) U uw, forms a closed PL curve, as Figure [5] shows. By 
Fenchel's theorem we have 

T K (L(wu) Ural) > 2tt. (7) 

Denote the exterior angle of the PL curve L(wu) U v/w at w as cti (Figure 
that is, 

«i = r](ml}, wv\). 

By Inequality [4], 

Since mb C fi, we have that rj(ui!j,nn) = |. Note also that ^(n^, MWi) = 
0(t o ). So 

«1 < | + 0(*o)- 

Denote the exterior angle of the PL curve L(wu) U » at w as a<i- By the 
definition of exterior angles, we have a 2 < 7T, so that 

T K (L(wu) Uri)=«i + T K (L(W)) + a 2 < ^ + 0(*o) + T K (L(W)) + vr. 

It follows from Inequality [7] that 

^ + fl(t )+T K (L(wm)) + 7r>27r, 
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so 



T K (L(wu))+9(t ) > |. (8) 



By L{wu) C L r , we have 



7T 



T K (L r ) + 9(to) > T K (L(wu)) + 9(to) > 
But this contradicts Inequality [6j □ 

Lemma 4.4. Any plane parallel to fl intersects L at no more than a single 
point. 

Proof: Suppose fl is a plane parallel to fl. If fl fl L = 0, then we are done, 



so we assume that fl fl L ^ 0. If fl = fl, then Lemma 4^3 applies, so we also 
assume that fl ^ fl, implying that w ^ fl. 

Consider two closed half-spaces EI; and EL such that EI; U M r = K 3 and 
Mi n H r = fl. Since fir\Li = fir\L r = w and L = L; U L r is simple, we can 
assume that Li C EI; and L r C i r 

Suppose without loss of generality that fl C il r , as shown in Figure [6] 
Then since L\ C EI/ and H| n H r = fl ^ O, we have fl C\ Li = 0. Since we 
assumed fl fl L ^ 0, it follows that fl L r 7^ 0. Now, it suffices to show that 
fl H L r is a single point. 

Since L r is compact and oriented, let w denote the first point of L r , at 
which fl intersects L r . Since fl \\ fl and f2 7^ f2, we have w ^ w. We shall 
show that fid L r = w. 

Denote the sub-curve of L r from its initial point vq to w as K\, and the 
sub-curve from w to its end point v n as i^2, as shown in Figure [6] Since w is 
the first intersection point of fl fl L r , but K\ ends in u>, then it is clear that 
fl fl K\ contains only w. Then in order to show fl fl L r = w, it suffices to 
show that fl fl K 2 = w. 

If w = v n , then it is the degenerate case: K 2 = w, and we are done. 
Otherwise, there is a vertex v k for some k G {1, . . . , n} such that is the 
non-degenerate initial segment of K 2 , where w 7^ v k . Now we shall establish 
the inequality: 

— — y 7r 

T K {K 2 ) +r](n n ,wv k ) < -, 

to guarantee a single point of intersection, similar to arguments previously 
given in Lemma 4.3| To this end, we use Inequality [4] to note that 



r)(n n , wv k ) < i](n n , vtfv{) + r](v^, wv k ) = 6(t ) + r](v^, wv k ). (9) 
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L= Li U Lr Lr=KiUK 2 




Figure 6: A parallel plane intersecting L 



The proof will be completed if we can show that 



T K (K 2 ) + 9(t ) + r](vtfVi, wv k ) < 



71 



(10) 



Casel: The intersection w is not a vertex, that is, w ^ v k -\. Then 
w is an int erior jp oint of v k -\V k , and , hence T K (Ki) = rj(v ovl, viv 2 ) + • • • + 



r}(v k - 2 Vk-i,Vk-iw), and r](v k _iw,wv k ) = 0. By Lemma 



4.2 



r}(vtfvi,wv' k ) < r](v^,v 1 ~v 2 )+. . .+r)(v k - 2 v k -i,v k -iw)+r)(v k - 1 w, wv k ) = T K {K X 
So 



7» „"T-^ 



~± — 



T K (K 2 ) + 9(t ) + V (^Ul,wv k ) < T K (K 2 ) + 9(t ) + T K {K X ). 
We also have 



-± - 



T K (L r ) = T K {K X ) + n(v k - x w,wv k ) + T K (K 2 ) = T K (K X ) + T K (K 2 ), 



(since r](v k -iW, wv k ) =0), so that 

T K {K 2 ) + e{t ) + ritwti, ^Z) < T K (L r ) + 9(t ) 
which is less than |, by Inequality M 
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Case2: The intersection w is a vertex, that is, w = Vk-i, then T K {K\, 
r}(vtfvi, v^V2) + • • • + i](v k _ 3 v k _ 2 , v k - 2 w). By Lemma 



4.2 



< T](vtfvi, v^vt) + • • • + ri(v k -3Vk-2, v k - 2 w) + r](v k - 2 w, wv k ) 

< T K {K\) + i](v k ^ 2 w, wv k ). 



So 



T K (K 2 ) + 0(i o ) + rj(v vi, wv k ) < T K (K 2 ) + 0(t ) + T R (K\) + r](v k _ 2 w, wv k ). 
But by the definition of the total curvature for a PL curve, 



i — 



^(#2) + + V(Vk-2W, wv k ) = T K (L 



So 



T K (K 2 ) + 6(t ) + 7i(wft, wv k ) < T K (L r ) + 6(t ) 
which is less than |, by Inequality [6] 



So Inequality [TO] holds, whi ch is an inequality analogous to Inequality |6j 
If in the proof of Lemma 4.3, we change Q t o Q, L r to K 2 and 9(to) to 
r](nQ,wv k ), then a similar proof of Lemma 4.3 will show that Q n = 
This completes the proof. □ 



Lemma 4.5. For an arbitrary t G [0,1], i/ie cfoc D r (t Q ) intersects C at a 
unique point, and also intersects L at a unique point. 



Proof: First, we have C(to) G D r (to) fl C. If there is an additional point, 
say C(ti) G D r (to) fl C where t\ ^ to, then we have that C(ti) ^ C(to) 
because C is simple, and hence D{t\) ^ D(to). Since also C(t\) G D r (ti), 
we have that C(ti) G D r (t ) fl D r {ti). But this contradicts the non-self- 
intersection of T. So D r (t ) fl C must be a unique point. 

Now, we show that D r (t ) fl L 7^ 0. If £ = or t = 1, then since 
L(0) G D r (0) and L(l) G D r (l), we have that D r (t ) H L ^ 0. 

Otherwise if to £ (0, 1), then assume to the contrary that D r (to) fl L — 0. 
Since L C T by Hypothesis 1, the contrary assumption implies that L C 
T \ D r (t ). Because C is an open curve, we have that D r (0) 7^ D r (l). So 
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r \ D r (t ) consists of two disconnected components, but this implies that L 
is disconnected, which is a contradiction. So 



D r (t )nL^ 



(11) 



Since D r (to) || f2 (as discussed in Section 4.1 ), Lemma 4.4 implies that the 



plane containing D r {to) intersects L at no more than a single point, which, 
of course, further implies that D r (t ) intersects L at no more than a single 
point. This plus Inequality 11 shows that D r (t ) D L is a single point. 

If -D r (t ) H L = D r (ti) n L for some £i 7^ £o> then D r (t ) and D T (t{) 
intersect, which contradicts the non-self-intersection of V. So there is an 
one-to-one correspondence between the parameter t and the point D r (t) n L 
for t G [0, 1], which shows the uniqueness. □ 



Lemma 4.6. The map L(t) given by Equation\Qis well defined, one-to-one 
and onto. 



Proof: It is well defined by Lemma |4.5| Suppose L(t\) = £(£2), then 



D r (ti) fl L = D r (t 2 ) fl L which is not empty by Lemma 4.5 So D r {t\) fl 
D r (t 2 ) 7^ 0. Since T is nonsingular, it follows that D r (ti) = D r (t 2 ). Since 
C is simple, if D r {ti) = D r (t 2 ), then t\ = t 2 . Thus L is one-to-one. Since 
L C T, each point of L is contained in some disc D r (t). So L is onto. □ 

Lemma 4.7. The map L{t) given by Equation is continuous. 

Proof: Let r tlt2 be the portion^] of T corresponding to [£1,^2], that is 

r tlt2 = |J D r (t). 

te[ti,t 2 ] 



Suppose that s G [0, 1] is an arbitrary parameter. Then by Lemma 4.6 , there 
is a unique point q G L such that q = L(s) = D r (s) fl L. We shall prove the 
continuity of L{t) at s by the definition, that is, for Ve > 0, there exists a 
5 > such that \t — s\ < 5 implies \ \L(t) — L(s)\\ < e. 

Note that D r (s) divides T into r 0s and T s i. Since C is an open curve, 
it follows that D r (0) 7^ D r (l), and that r 0s and T s i intersect at only D r (s). 



By Lemma 4.5, D r (s) fl L is a single point, so L is divided by D r (s) into two 
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Figure 7: If \s — r| < S, then \\q — q' 



< € 



sub-curves, denoted as K\ and K 2 , that is K\ C T 0s and C T s i, as shown 
in Figure [7} 

Casel: The parameter s is such that s / and s ^ 1. Consider T s i 
first. Since K 2 is oriented, we can let v be the first vertex of K 2 that is 
nearest (in distance along K 2 ) to q. For any < e < \\qv\\, let q' & qv such 



that ||gg'|| = e. By Lemma 4.6 g' = £(t) = D r (r) fl L for some r G (s, 1]. 

First, note qq' fl intr sT 7^ 0. To verify this, observe qq' Cqv C K 2 C F s i 
and T s i = r sT U r r i, so qq' C T ST U r r i. If qq' fl intr sT = 0, then the segment 
qq' is contained in D r (s) U T rl which is disconnected. This implies qq' is 
disconnected, which is a contradiction. 

Secondly, note that the subset r sr of a nonsingular pipe section is con- 
nected (since C is C 1 ), and qq' is a line segment jointing the end discs of T ST , 
and has intersections with interior of r sr . This geometry implies that 

g? c T ST . (12) 



Let 5 = t — s. For an arbitrary t & (s, s + 8) = (s, r), Inclusion 12 implies 
that !/(/!;) = D r (t) fl gg'. Since neither L{t) ^ q or L(t) 7^ g', it follows that 
L(t) E int(gg'). So 

\\L(t) - L(s)\\ < ||g?|| =e. 

This shows the right-continuity. We similarly consider the T 0s and obtain 
the left-continuity. 



If t\ = t2, then r tl t 2 is a single disc. 
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Case2: The parameter s is such that s = or s = 1. We similarly 
obtain the right-continuity if s = 0, or the left-continuity if s — 1. □ 



Theorem 4.1. If L and C satisfy Hypotheses 1 and 2, then the map h 
defined by Equation\^is a homeomorphism. 



Proof: By Lemma 4.6 L(t) is one-to-one and onto. By Lemma 4.7 L{t) is 
continuous. Since L is defined on a compact domain, it is a homeomorphism. 

Note that C is simple and open, so C{t) is one-to-one, and it is obviously 
onto. The map C(t) is also continuous and defined on a compact domain, 
so C(t) a homeomorphism. Since h is a composition of C _1 and L, h is a 
homeomorphism. □ 



Remark 4.3. A very natural way to define a homeomorphism between sim- 
ple curves C and L would be by f(p) = L(C^ 1 (p)) . An easy method to extend 
f to a homotopy is the straight-line homotopy. However, we were not able to 
establish that a straight-line homotopy based upon f would also be an isotopy, 
where it would be necessary to show that each pair of line segments generated 
is disjoint. Our definition of h in Equation\^was strategically chosen so that 
this isotopy criterion is easily established, since the normal discs are already 
pairwise disjoint. 



5. Construction of ambient isotopies 



Note that L and C fit inside a nonsingular pipe section (Remark 3.1) 
T of C. For a similar problem, an explicit construction has appeared 



Section 4.4] [TTJ. The proof of Lemma 5.3, below, is a simpler version of 
a previous proof [HI Corollary 4]. The construction here relies upon some 
basic properties of convex sets, which are repeated here. For clarity, the 



complete proof of Lemma 5.3 is given here 



Lemma 5.1. 177], Lemma 6] Let A be a compact convex subset of M 2 with 
non-empty interior. For each point p G int(A) and b G OA, the ray going 



from p to b only intersects dA at b (SeenFigure 8(a).) 



6 The Images in Figure [8] were created by L. E. Miller [18] and are used, here, with 
permission. 
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(a) Rays outward 



(b) Variant of a push 



Figure 8: A convex set and a map defined on it 



Lemma 5.2. 177], Lemma 7] Let A be a compact convex subset of M 2 with 
non-empty interior and fix p G int(A). For each boundary point b G dA, 
denote by \p,b] the line segment from p to b. Then A = U&eaA[P>^]- 

Lemma 5.3. There is an ambient isotropy between L and C with compact 
support ofT, leaving dT fixed. 

Proof: We consider each normal disc D r (t) for t G [0, 1]. Let p = D r (t)DC 
and q = h(p) with h defined by Equation [3j then define a map F Piq : D r {t) — > 
D r (t) such that it sends each line segment \p,b] for b G dD r (t), linearly 



onto the line segment [q, b] as Figure 8(b) shows. The two previous lemmas 



(Lemma 5.1 and Lemma 5.2 ), will yield that F p>q is a homeomorphism, leaving 
dD r {t) fixed [W\ Lemma 4.4.6]. 

In order to extend F p ^ q to an ambient isotopy, define H : D r (t) x [0, 1] — > 
D r (t) HH Corollary 4.4.7] by 

H(v,s) = 



(1 — s)p + sq if v = p 

•^*p,(l— s)p+sq \P 



where F Pj ( 1 _ s ) p+S q is a map on D r {t) analogous to F Ptq , sending each line 
segment [p, b] for b G dD r (t), linearly onto the line segment [(1 — s)p + sq, b]. 

It is a routine [HI Corollary 4.4.7] to verify that H(v,s) is well defined 
on the compact set D r (t), continuous, one-to-one and onto, leaving dD r (t) 
fixed. Now, we naturally define an ambient isotopy T t : 1R 2 x [0, 1] -> M 2 on 
the plane containing D r (t) by 



T t (v, s) 



H{v,s) if veD r (t) 
v otherwise. 
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We then define T : R 3 x [0, 1] ->■ M 3 by 



T(v,s 



T t (v,s) if veD r (t) 
v otherwise. 



The fact that the normal discs D r (t) are disjoint ensures that T is an ambient 
isotopy [TBI Corollary 4.4.8], with compact support of T, leaving dT fixed. □ 



Now we apply this result to a simple, regular, composite, C 1 Bezier curve 
B and the control polygon V . As discussed previously, there exists a non- 
singular pipe surface [17] of radius r for B, denoted as S#(r). Denote the 
nonsingular pipe section determined by 5g(r) as IV Also, for each sub- 
control polygon of B, there exists a corresponding nonsingular pipe sections. 
Denote the nonsingular pipe section corresponding to the fcth control polygon 

as r fe . 

Theorem 5.1. Each sub-control polygon P k of a Bezier curve B will even- 
tually satisfy Hypotheses 1 and 2 via subdivision, and consequently, there will 
be an ambient isotropy between B and V with compact support ofT&, leaving 
dTjs fixed. 



Proof: By Lemma 4.1, Hypotheses 1 and 2 can be achieved by subdi- 
visions. Now consider each sub-control polygon P k satisfying Hypotheses 1 



and 2, and the corresponding Bezier sub-curves B k . Use Lemma 5J3 to de- 
fine an ambient isotopy : IR 3 x [0, 1] — > 1R 3 between B k and P k , for each 
k E {1,2,..., 2*}. Define * : 1R 3 x [0, 1] -> 1R 3 by the composition 

$ = $lO$ 2 o...o \]> 2l . 

Note that fixes the complement of int(r&), and int(rfe) fl int(Tfc/) = 
for all k k'. So the composition \1/ is well defined. Since each ^/^ is an 
ambient isotopy, the composition \1> is an ambient isotopy between B and V 
with compact support of Tg, leaving 3Yq fixed. □ 

6. Sufficient subdivision iterations for ambient isotopies 

Now we consider sufficient numbers of subdivision iterations to achieve 



the ambient isotopy defined by Theorem 5.1 By Theorem 5.1, we shall 



have a control polygon that satisfies Hypotheses 1 and 2. The number of 
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subdivisions for Hypothesis 1 is given in the companion paper [TJ)1 Lemma 
6.3]. To obtain the number of subdivisions for Hypothesis 2, we consider the 
following lemmas, for which we let V'(t) = l'(P,i)(t) (the first derivative of 
the control polygon V), and denote the angle between B'{t) and V'(t) as 9(t), 
forte [0,1]. 

Lemma 6.1. For any < v < | ; let N(u) be the value given in the 
companion paper [15], Equation 20], then after N(v) subdivisions, we have 
max te[0)1] 6(t) < v. 

Proof: Note (similarly as [H2, Inequality 6]) 



\\B\t)\\.\ 


B'(t) ■ 1 
l-co sW )) = l- ||B , (( ^ 

\V'(t)\\-V'(t)-V'(t)+V'(t)- 


o'(t) 
\V'(t)\\ 

V'(t) - B'{t) ■ V'(t) 




\B'(t)\\-\ 


\V'(t)\\ 



\mt)\\-\ 


\v(t)\\ 




mm 



< ■■- \\B'(t)-V'(t)\\ 2\\B'(t)-V'(t)\\ 

\\B'{t)\\ ~ 

where a = min{||i3'(t)|| : t G [0,1]} (Recall a > 0.) From the companion 
paper [T5J Inequality 12] 

max\\B'(t)-V'(t)\\<B> dist (i), 
te[o,i] 

we have 

l-cos(0(t)) < 2B ' d ^ l \ 
a 

Then from the companion paper [T51 Inequality 16] (Note that this inequality 
is valid for i > N(u)) 

2(2^(0+ 7 /(ny)) _ 

--BUN.) Kl C ° S{Uh 

we get 

i - cos Wt)) < < jgS4±^) < , _ cosM . 

So cos(#(£)) > cos(u). Note that 8(t) is an exterior angle with measure be- 
tween and it and < v < | by the hypothesis. Since cos a; is monotonically 
decreasing when x G [0, |], we have 8(t) < v. This is true for all t G [0, 1], 
so maxt e [o,i] 9(t) < v. □ 
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Lemma 6.2. Performing N(j-) subdivisions will yield a control polygon V 
with each sub- control P k satisfying Hypothesis 2: 

T K (P k ) + max 6(t) < J. 

Proof: Recall that N(j-) given in the companion paper [T5l Equation 
20] is a sufficient number of subdivisions to produce V such that all exterior 
angles are less than j-, as shown in the companion paper [Tol Theorem 6.1]. 
Note that each P k has n — 1 exterior angles. So after N(-^) subdivisions, we 
have 

By Lemma 6.1, after iV(^) subdivisions, we also have 

71 

max 9{t) < — . 

te[o,i] 2n 

So, after N(j-) subdivisions, we obtain 

T K {P ) + max t < -\ ^ + — = - . 

te[o,i] 2n 2n 2 

□ 

Let 

iV* = max{iV(|-),jV'(r)}, (13) 
where r is the radius of S r (B). 

Theorem 6.1. Performing N* or more subdivisions, where N* is given by 



Equation 13, will produce an ambient isotopic V for B. 



Proof: According to the companion paper [Td"| Lemma 6.3], Hypothesis 1 



is satisfied after N'(r) subdivisions. By Lemma 6.2, Hypothesis 2 is satisfied 
after N{^-) subdivisions. Then Theorem 5.1 can be applied to draw the 



conclusion. □ 
Hypothesis 2 is motivated by the weaker condition T K (P k ) < |. We 
couldn't derive the same results by using this weaker condition instead, but 
our Hypothesis 2 requires at most one more subdivision, as shown below. 
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Remark 6.1. The sufficient number of subdivisions for Hypothesis 2 is at 
most one more than that for the weaker condition T K (P k ) < |. 



Proof: By the above proof of Lemma 6.2 we know that iV( 2 ,^ ) subdi- 
visions are sufficient to guarantee T K (P k ) < |, while N(j-) subdivisions are 
sufficient to obtain T K (P k ) + max te p ^ 9(t) < |. The function N(-) is given 
in the companion paper [T5J Equation 20], that is, 

N(v)=max{N 1 ,N 2 ,]og(f(v))}, 

where Ni and N 2 are constants for a given B, and the function /(■) is defined 
in the companion paper [T5J Equation 19]: 

fW) = i , 

1 — COS V 

where eg is a constant. Note that, 
and 

iV(|-) = max{iV 1 ,iV 2 ,log(/(^))}. 

So comparing N( 2 ,*_^ ) with N(^) is essentially comparing log (/( 2(^1) )) 
with log(/(^)). So we consider 

7T„ , * „ , /(£) _ l-COs(^) 



l0g(/(— )) - log(/(- -)) = log -j^^ = log 



*2n" owv 2(n-l)" /(^) 1 - cos(^) ' 

As before, we consider only n > 3, since for B of degree 1 or 2, the control 
polygon is trivially simple and ambient isotopic to B, provided the regularity 
ofB. 

Consider the expression inside the logarithm as a function of x: 

1 " C0S (2(^T)) 



i-cos(^; 



Since its derivative is negative for x > 3, we have that 

1 — cos 



■<*(/(£)) - log(/(^, )) < log ' I ;^ ( 'g > - - 0.78. 
This implies that 

N( — )-N( , 71 r ) < 1, 

v 2n ; v 2(n-l) ; 
for any n > 3. □ 
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7. Conclusion 

We showed that sufficiently many subdivisions produce a control polygon 
ambient isotopic to a given Bezier curve B. We established closed-form for- 
mulas to compute a priori sufficient number of subdivisions to achieve the 
isotopy. The proof here is based on explicitly constructing homeomorphisms 
and ambient isotopies, which is motivated by the goal of minimizing sub- 
division iterations for practical applications. In contrast, the pure existence 
proof [T2] requires the convex hulls of sub-control polygons to be contained in 
the tubular neighborhood determined by a pipe surface and may need more 
subdivision iterations and produce too many PL segments. This method 
here has another theoretical consequence: we established two criteria (Hy- 
potheses 1, 2) to ensure the ambient isotopy between a PL curve and a C 2 
curve connected at the end points. 

These results have potential applications in computer graphics, computer 
animation and scientific visualization, especially in visualizing molecular sim- 
ulations. 
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